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Abstract 

We perform a Monte Carlo Renormalization Group analysis of the critical be- 
havior of the ferromagnetic Ising model on a Sierpihski fractal with Hausdorff di- 
mension df ~ 1.8928. This method is shown to be relevant to the calculation of 
the critical temperature T c and the magnetic eigen-exponent yu on such structures. 
On the other hand, scaling corrections hinder the calculation of the temperature 
eigen-exponent yt- At last, the results are shown to be consistent with a finite size 
scaling analysis. 
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75.40.Mg Numerical simulation studies 
89. 75. Da Systems obeying scaling laws 

1 Introduction 

Since Mandelbrot attracted people's attention to the fractals in the 60's and 70's pQ, 
scientists began to take auto-similarity and scaling invariance as fundamental rules of 
nature. In the description of second order phase transitions, Widom's homogeneity hy- 
pothesis |2] and Kadanoff 's real space renormalization group scheme [3] are based on the 
invariance of the physical behavior under any change of scale at the criticality where the 
correlation length is divergent. The validity of these hypothesis and scheme has been 
verified in different systems with an integer dimension which possess a translational sym- 
metry As a matter of fact, a system with translational symmetry is auto-similar and 
can be considered as a particular case of fractal. Fractals are natural candidates to repre- 
sent systems with non-integer dimensions. It turns out to be of fundamental relevance to 
know if the hypothesis and the scheme mentionned above still work on a general fractal 
system, where the translational symmetry is lost and replaced by scale invariance. 

The critical behavior of the Ising model on fractal lattices has been firstly studied 
by Gefen and his co-workers in the 80's Their works were based on the Migdal- 
Kadanoff bond-moving renormalization method and showed that the topological features 
of the fractals plays an important role in the determination of the critical behavior. As a 
main result, they found that a phase transition at non-zero temperature can occur only 
if the fractal has an infinite ramification order. Later, Bonnier et al. used an alternative 
decimation method P and high-temperature expansions [7j to study the phase transition 
on various fractals, namely Sierpihski carpets. They also found that the critical exponents 
depend on the geometrical properties of the fractal. In the above studies, the applied 
methods are approximative and the results are not always consistent. Moreover, in these 
theoretical analyses, the spins were placed at the corners of the occupied squares of 
the Sierpihski carpets; in these cases, the number of spins doesn't follow a power law 
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of the lattice size and the Hausdorff dimension should not be expected to enter in the 
description of the critical phenomena on such systems. Recently, due to the progress in 
the simulation methods and the growth of the computer power, the critical behavior of 
the Ising model on fractals of Hausdorff dimension df between 1 and 3 has been studied 
numerically in a much more precise way [HI HH IUJ1 E] • The results showed that the finite- 
size scaling (FSS) analysis works in the case of fractals, although the convergence towards 
the thermodynamical limit can be very slow when df < 2, and that the hyperscaling law 
df = 2/3 /u + 7/V is verified. Moreover, discrepancies with the predictions of the e- 
expansions ^3] were observed. The universality of phase transitions on fractals is said 
to be weak [5J [J3]. It is worth noting that the spins were placed at the center of the 
occupied squares of the Sierpihski carpets in these studies; consequently, the number 
of the spins increases as a power law of the lattice size with an exponent equal to the 
Hausdorff dimension df. In this case, df takes the place of the space dimension in Widom's 
homogeneity hypothesis because of the conservation of the total free energy under a 
change of length scale; the generalization of this hypothesis in the case of fractals is quite 
straightforward [T2"] . 

Monte Carlo renormalization group (MCRG) has been shown to be a powerful tool 
in the study of critical phenomena [To"! ITdI IT71 ITS] ITU] . It combines Monte Carlo (MC) 
simulation techniques with the analysis of the real space renormalization group (RSRG). 
It has been, so far, applied to systems with translational symmetry. A RSRG analysis in- 
volves an infinite number of couplings when dealing with the Hamiltonian. Theoreticians 
usually truncate the number of couplings to a finite one in order to make the calculations 
tractable. It, therefore, makes RSRG an approximate method and the accuracy of the 
results is linked to the number of the couplings considered. Even if only few couplings 
are considered, the calculations associated with the process of decimation remain tedious. 
With the aid of MC simulation, a renormalization-group calculation with a larger number 
of couplings becomes possible and the results are expected to be more reliable. However, 
some crucial problems are encountered in a MCRG study: 

(1) The critical slowing down: 

Hypothetically, the accuracy of the thermodynamical averages can be improved by 

3 



performing a larger number of MC steps to reduce the statistical errors. In fact, due 
to the divergence of the correlation length at the criticality, local MC algorithms 
(for e.g., the Metropolis one) suffer from critical slowing down, which hinders the 
computation of accurate thermodynamical averages. Fortunately, the use of cluster 
algorithms, for e.g., the Swendsen-Wang algorithm j20] and the Wolff algorithm 
[21], enables to improve significantly the efficiency of the simulations. 

(2) The finite size effects: 

RSRG method involves an infinite size lattice, whereas computer simulations can 
only be performed on finite size lattices. One should perform the simulation on 
several large systems and then extrapolate the results to the infinite one. 

The MCRG method, according to which one plays the game of changing the length scale, 
should be a natural approach to explore the critical behavior of phase transitions on 
hierarchical lattices. The purpose of this paper is to study the critical behavior of the 
Ising model on a Sierpihski carpet using the MCRG method. As far as we know, it is 
the first time that the MCRG method is used in the case of fractal lattices. The paper 
is organized as follows : In Sec. II, we present the model and recall briefly the MCRG 
method. We explain how the simulation is setup in Sec. III. Sec. IV is devoted to the 
calculation of the critical temperature T c , the temperature eigen-exponent y t , and the 
magnetic eigen-exponent yt- A FSS analysis is set out in Sec. V for a consistency check. 

2 Sierpiriski carpet and MCRG method 

A 3 by 3 square lattice with its central subsquare removed is chosen as "generating cell" . 
The Sierpihski carpet we deal with is constructed iteratively from this generating cell: 

1. We take the generating cell as the carpet at the first iteration step and denote it 
SC(3,1,1). 

2. The carpet at the (A; + l)-th iteration step SC(3, 1, k + 1) is constructed by enlarging 
three times the size of the carpet at the k-th iteration step and replacing each 
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enlarged un-removed subsquare by the whole generating cell. This carpet becomes 
a "true" fractal in the mathematical sense when k tends to infinity. 

The size of the lattice SC(3, 1, k) is equal to = 3 k . A spin is set at each center of an 
un-removed subsquare (so-called occupied square). The number of spins on SC(3, l,k) is 
therefore equal to Nk = 8 k . We notice that is equal to Lk df where df = log 8/ log 3 ~ 
1.892789 is the Hausdorff dimension. The critical behavior of the ferromagnetic Ising 
model on this Sierpihski carpet has been recently studied independently by P. Monceau 
et al. [HIE] and J. Carmona et al. jH] by means of MC simulation. Thus, estimations of 
the critical temperature are available and the comparison between the two methods will 
be possible. 

Let us briefly recall the MCRG method. The geometrical structure of a Sierpihski 
carpet remains invariant during a RSRG process involving a change in the length scale 
from 1 to b, provided that b is equal to an integer power of the size of its generating 
cell. In our case, we set b equal to 3 and increase the length scale by this factor at each 
renormalization step. The majority rule is used to decimate the spin blocks: a new spin 
is assigned to a given block according to the sign of the summation of the spin states of 
the 8 occupied sites in this block; if the sum is zero, a spin state +1 or —1 is assigned to 
the block with the same probability. The reduced Hamiltonian after n renormalization 
steps {n is a positive integer) reads as 



ft (n) ({^i n) }a=i) = E K[ : )s{ : ] 

a=l 



where {i^i }^=i ^ s the coupling set after n renormalization steps and S^'s are the 
conjugate lattice sums of the spin products on the reduced network. The matrix T, which 
describes the flow of the couplings from the n-th renormalization step to the (n + l)-th 
one, is defined by 



and satisfies the following relation: 

oo 

[os}% +1 ' n) = £ [ds}% +1 ' n+ v r^ n \ (3) 

Q = l 



^(n+1) 
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where [dS] { ^ m) = d(S^)/dK^ ] can be calculated by (S^S^) - (S^)^). One 
can hence calculate the matrix T by inverting the above relation. Two critical exponents 
enable to describe the static scaling properties of the system: the temperature eigen- 
exponent y t and the magnetic one They describe respectively the scaling behavior 
of the reduced temperature t = T/T c — 1 (where T c is the critical temperature) and the 
external magnetic field h near the criticality under a change of the length unit from 1 to 
b: 

t t' = b yt t , 

h -> h' = b Vh h . 

b Vt and b Vh are associated with the two relevant directions of the renormalization flows 
and can be obtained by finding the largest eigenvalues of the T matrix in the even- and 
odd-coupling subspaces, respectively. 

3 Simulation setup 

Up to 4-spin couplings have been considered in our study and the interaction range has 
been restricted within a 3 x 3 square. There are 25 even couplings and 11 odd couplings as 
shown in Fig. 1(a) and Fig. 1(b), respectively. The symmetry number obtained by rotating 
and reflecting a given coupling is indicated in the third column of the figure. These 
couplings are listed in decreasing order according to the importance factor proposed by 
Blote et al. [17] : F = (2 a / 2 f)~ 1 , where r is the average distance between the s spins. 
The value of F is given in the right column of the figure. The simulations have been 
carried out at two different temperatures, T = 1.4813 obtained by Carmona et al. jHj and 
T = 1.4795 obtained by Monceau et al. [II] . The size of the starting Sierpihski carpet 
varies between 3 4 = 81 and 3 8 = 6561. The simulation procedure is organized as follows: 

1. An equilibrium spin configuration on the starting Sierpihski carpet SC(3, l,k) is 
generated by the Wolff algorithm [21] at the simulation temperature. 

2. The lattice sums conjugated to the 25 even couplings and to the 11 odd ones are 
calculated. 
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3. A decimation is done, according to the majority rule, by dividing the size of the 
Sierpihski carpet by 3 and the lattice sums associated with the renormalized carpet 
are calculated. 

4. The step 3 is repeated until the size of the renormalized Sierpihski carpet is smaller 
than 9. 

A data sample of the lattice sums is built up from 10 5 steps of the above simulation, 
with periodic boundary conditions. At each simulation temperature, for the five different 
sizes of the starting Sierpihski carpet, 10 independent data samples are collected. The 
errors bars are estimated from a standard statistical analysis of these data. 



4 MCRG Results 

4.1 Critical temperature T c 

The size of the lattice after n + p renormalization steps from the starting Sierpihski 
carpet of size Lk+ P is equal to the size of the lattice after n renormalization steps from 
the starting Sierpihski carpet of size At the fixed point K*, the thermodynamical 
averages of the lattice sums calculated from renormalized carpets with the same size 
should be independent of the size of the starting Sierpihski carpet. Let us assume that 
the initial coupling vector K^ ) is close to K*. The difference of the lattice sums on the 
reduced lattices of equal size, derived from the small displacement <5K = K* — ~KS°\ 
should satisfy the following relation: 

- (sir'ML - - £ (K + J!r 0) - [«&JS*) K( „, <«) 



where the suffix Lk+ P or Lk indicates that the physical quantity is obtained from the 
starting Sierpihski carpet of size Lk+ P or Lk- In the even coupling subspace, the nearest 
neighborhood coupling K\ represents the inverse of the temperature. If we neglect the 
contribution from the directions other than the nearest neighborhood coupling and set 
a = 1 in Eq.(j3J, the critical temperature can be calculated by 

Tf +P > k \n) = (KlY 1 =\Kf ] +5Ki' ] 



where the symbol T\ k+P ^ (n) indicates that the temperature is obtained from two starting 
carpets of size 3 k+p and 3 k after n + p and n renormalization steps respectively. Fig. 2 
shows the evolution of xJ fc +P' fc ] ( n ) calculated from pairs of the starting Sierpihski carpets 



at the simulation temperatures T sim = = 1.4813 and T sim = K\ 



-i 



,(0) 



-i 



= 1.4795. 



These results deserve the following comments : 

i) At the two simulation temperatures, T^ k+P ' k \n) decreases as n increases and tends 
to converge to some value. 

ii) If we fix p — 1, the value of T^ +P,k \n) decreases as k increases. We can see that 
Tf^(n) > Tf'%n) > T] 7 ' 6 l(n) > Tf^(n) and the curves T^^(n) and Tf^(n) are very 
close to each other. 

iii) Tf +P > k \n) decreases as k decreases if the size of the largest starting Sierpihski 
carpet 3 k+p is fixed. One can observe that Tf> 7 \n) > T^(n) > T^(n) > Tf A \n) 
where the size of the largest starting carpet is fixed at 6561. These curves tend to 
converge to some temperature while n increases. 

The best estimation for the critical temperature can be obtained from the 6-data 
points T] 8 ' 7 ] (n) for n — 0, 1, • • • , 5, by performing the following three-parameter (T c , A , n ) 
fit: 

T [k+P ' k \n) =T C + A 3- n/no . (6) 
The results of the fits at the two simulation temperatures lead to : 

at T sim = 1.4813, (T c , A , n ) = (1.47959(15), 0.02718(28), 0.839(25)) 

with the reliability R 2 = 0.99968 and 
at T sim = 1.4795, (T c , A , n ) = (1.47933(27), 0.02827(49), 0.868(42)) 

with the reliability R 2 = 0.99911 . 

The best value of the critical temperature can be estimated by taking the average of the 
results of the two fits; we hence have T c = 1.47946(16). 
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In the following subsections, the eigen-exponents will be calculated from simulations 
carried out at T sim = 1.4795, since it lies closer to the best estimate of the critical 
temperature T c = 1.47946(16). The simulations carried out at T sim = 1.4813 will provide 
a consistency check. It is well known that the exponents provided by the MCRG method 
are not very sensitive to small deviations from the critical temperature. 

4.2 Temperature eigen- exponent y t 

We work in the even-coupling subspace and consider in turn the first 1 to 25 couplings 
shown in Fig. 1(a). We calculate the T matrix in the subspace at each renormalization 
level by inverting Eq.fjSJ) with the help of the Gauss- Jordan elimination method (§2.1 in 
Ref.j221)- The dimension of T is N coupHng x N coupling where N coupHng varies between 1 
and 25. The largest eigenvalue At of T is obtained by applying the eigenvalue searching 
method for a real non-symmetric square matrix given in §11.5 and §11.6 of Ref.j22]- The 
temperature eigen-exponent is then calculated by y t = log 3 X t . Fig. 3 shows the evolution 
of the eigen-exponent y t as a function of the number of even couplings N coup i ing , obtained 
from starting Sierpihski carpets of different sizes and at different renormalization levels. 
One can see that y t tends to converge towards stable values when N coup n ng increases, 
within the statistical errors, excepted at the highest renormalization levels where the 
sizes of the reduced lattices are 3 3 and 3 2 . We estimate yt in the infinite couplings 
limit by taking the average of the values of yt obtained from 16 to 25 even couplings. 
The average value at each renormalization level is reported in Table where we have 
disregarded the unreliable results at the highest renormalization levels. We find that the 
finite-size effect is not significant because at the same renormalization level, the values of 
yt obtained from the starting Sierpihski carpets of different sizes are about the same. A 
similar situation has been observed in the study of the 3D Ising model on regular lattices 
fTJEI]. In order to extrapolate the value of yt on the starting Sierpihski carpet of infinite 
size, we plot in the Fig. 4 the value of Xt obtained at level n — {n + 1) from the starting 
carpet SC(3, l,n + 4) with respect to n. We find that, up to the renormalization levels 
which were performed, At seems not to be convergent. This behavior is quite different 
from the one observed in the case of the 3D Ising model on regular lattices fTJ ^] where 
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At tends to converge to some value while n increases. We, hence, cannot estimate the 
convergence by performing the three-parameter fit (A|,a t , u t ) proposed in Ref.|18j: 

A t = A* + a t 3~^ n . (7) 

The lack of convergence in the case of the Sierpihski carpet may be due to a small value of 
the scaling correction exponent Ut . At the present stage, we are only able to provide an 
upper bound for the Sierpihski carpet of infinite size: At < 1.781(10) or y t < 0.5254(51). 
The same situation is also observed at T sim = 1.4813. The difficulty in estimating the 
value of y t has already been observed in the study based on FSS analysis [TT]. According 
to Widom's homogeneity hypothesis, y t is equal to the inverse of the correlation-length 
exponent v. Monceau et al. ^1] have shown that the maxima values of the logarithmic 
derivatives $j (for i = 1,2), defined by 

where E is the total energy and M the total magnetic moment, are affected by scaling 
corrections. Instead of scaling properly as a power law, the maxima values of $j (for 
i = 1,2) appear to exhibit a slight concavity in a log- log plot with respect to the lattice 
size. Monceau et al. could only provide an upper bound for y t from their MC study: 
\jv < 1/1.565 ~ 0.639. Our estimation is consistent with their result. 

4.3 Magnetic eigen-exponent 

We now focus on the odd-coupling subspace. We consider in turn the first N coupHng 
odd couplings shown in Fig.l(b) where N coup i ing = 1, 2, . . . , 11. The T matrix in the odd- 
coupling subspace is obtained from Eq.fjSJ) by considering only the lattice sums conjugated 
to the odd couplings. We denote A^ the largest eigenvalue of this matrix, and the magnetic 
eigen-exponent is calculated by yh = log 3 A^. Fig. 5 shows the evolution of yh as a funtion 
of the number of odd couplings at each renormalization level obtained from the starting 
Sierpihski carpet SC(3, 1, k) with k = 4, 5, • • • , 8. We can see that yh tends to be stable 
as N coup i ing increases. Larger fluctuations are observed at the highest renormalization 
level where yh is extracted from the reduced carpets of sizes 3 3 and 3 2 ; we will disregard 
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these values. We estimate yh in the infinite couplings limit by taking the average value 
of them when 7 to 11 odd couplings are considered. The results are reported in Tabled 
The size effect of the starting Sierpihski carpet on the evolution of with respect to the 
renormalization level is more significant than it is on the evolution of y t . However, the 
value of yh at level n — (n + 1) from the starting carpet SC(3, 1, n + 4) tends to converge 
to some value as n increases (see Fig. 6). We hence perform the three-parameter fit 
suggested by Baillie et al. [TSj : 

X h = X* h + a h 3~^ n . (9) 

The fit yields (X* h , a h , uo h ) = (7.38298(66), -0.0737(12), 1.86(12)) with a reliability equal 
to R 2 = 0.99949. According to these results, the magnetic eigen-exponent yh takes the 
value 1.81973(9) and the associated scaling correction exponent u>h the value 1.86(12). 
yh is slightly larger than the one obtained by Monceau et al. [TT]. They found that 
j/u — 1.732(4) and (3/v — 0.075(10). According to the relations: j/u = 2y h — df 
and j3/v = df — y h [I2] 5 their results correspond to yh = 1.812(2) and yh = 1.818(10), 
respectively. No significant discrepancy can be brought out from these results, since the 
relative difference between the exponent provided respectively by the MCRG method and 
MC simulations remains smaller than 0.5 percent. 



5 Results from FSS analysis 

The lattice sums associated with the first odd and the first even couplings represent 
respectively the total magnetic moment M and minus the total energy E. We can, 
hence, perform a FSS analysis for a consistency check. According to this analysis, the 
thermodynamical average of the total magnetic moment at T c should scale as (|M|) ~ 
L d f -/3/u = L y h and the Privative of the Binder's cumulant U = 1 - (M 4 )/(3(M 2 ) 2 ) at 
the critical temperature T c scales as 

^- = -(1 - 10 ($ 4 " 2$ 2) ~ L l '» = L yt , (10) 
a pb 

where L is the lattice size and (3b = {ksT)~ l is the inverse of temperature. In Fig. 7, we 
plot the values of (\M\) and dU/d 0b at T — 1.4795 for the 5 different lattice sizes from 
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3 4 to 3 8 in logarithmic coordinates. We find that the values of (\M\) line up along straight 
lines, with the lattice size covering several order of magnitude. The slope provided by a 
least-square fit from the 5 above points, is 1.8198(11), with a reliability R 2 = 0.99994. It 
turns out that this value is consistent with the value of yh obtained from MCRG method 
in the subsection 14.31 The behavior of dll/d^B is not affected by scaling corrections, 
(see Fig. 7), although the $j (for i = 1,2) exhibit a slight concavity as a function of 
the lattice size in a log-log plot. The reason why dU jd^s is not affected by scaling 
corrections may be a kind of "magic" cancellation of these scaling correction effects in 
the difference between $ 4 and 2$ 2 . The slope measured from dU/dfls is 0.449(6) with 
the fitting reliability R 2 = 0.99948. This value is consistent with the result obtained in 
the subsection 14.21 where y t < 0.5254(51). 

6 Conclusion 

The MCRG method has been shown to provide reliable values of the critical temperature 
T c and the magnetic eigen-exponent yh, for the Ising model in a case where the underlying 
network has a fractal structure. The scaling correction exponent uoh associated with 
yh can also be calculated. Moreover, difficulties are encountered in the estimation of 
the temperature eigen-exponent y t where, unlike regular lattices, the convergence of the 
renormalization flow in this direction seems to be very slow; nevertheless, an upper bound 
for y t can be provided. A similar slowness in the convergence of the logarithmic derivatives 
of the magnetization $j to the thermodynamical limit has already been brought out from 
MC simulations in the case of the same Sierpihski fractal jU] [TT]. In reference 11, scaling 
corrections have been shown to affect strongly the behavior of the $j's and only an 
upper bound for y t has been calculated from their finite size behavior. The difficulties 
in calculating y t for fractal structures arise in the two methods, but emerge differently. 
At last, the results provided by the MCRG method have been shown to be consistent 
with the one obtained by Monceau et al. ^T] and with an additional FSS analysis of the 
MCRG results. 
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Figure captions 

Fig.l (a) 25 even couplings and (b) 11 odd couplings, considered in our MCRG study. 
They are ordered according to the Blote's importance factor F whose value is given 
in the right colomn. The symmetry number is given in the next to last right colomn. 

Fig. 2 Evolution of Tj fc+P ' fc ](n) with respect to the different pairs of the starting Sierpihski 
carpets 5(7(3, 1, k+p) and 5(7(3, 1, k) (denoted by [k+p, k}) obtained at (a) T sim = 
1.4813, (b) T sim = 1.4795. 

Fig. 3 Variation of y t with respect to the number of even couplings N coup u ng at dif- 
ferent renormalization levels from the starting Sierpihski carpet 5(7(3, with 
= 4, 5, • • • , 8. 

Fig. 4 \ t obtained at level n — (n + 1) from the starting carpet 5(7(3, 1, n + 4). 

Fig. 5 Variation of with respect to the number of odd couplings N coup n ng at differ- 
ent renormalization levels from the starting Sierpihski carpet 5(7(3, 1, k) with k = 
4,5,- .-,8. 

Fig. 6 Xh obtained at level n — (n + 1) from the starting carpet 5(7(3, 1, n + 4). 

Fig.7 (|M|) and dU/d(3 B at T = 1.4795 on SC(3,l,k) with fc = 4, 5, • • • , 8 in logarithmic 
coordinates. 
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Table 1: Average value of y t calculated from the starting Sierpihski carpets at different 
renormalization levels. It is obtained by considering 16 to 25 even couplings at the 
simulation temperature T sim = 1.4795. At the highest renormalization level, these average 
value are disregarded (denoted by the symbol " — ") because the eigen-exponent is not 
stable with respect to the number of the even couplings. 



level 


SC(3,1,8) 


SC(3,1,7) 


SC(3,1,6) 


SC(3,1,5) 


SC(3,1,4) 


0-1 


0.7497(17) 


0.7437(20) 


0.7474(14) 


0.7506(26) 


0.7483(26) 


1-2 


0.7186(27) 


0.7235(14) 


0.7225(19) 


0.7221(17) 




2-3 


0.6465(26) 


0.6459(32) 


0.6466(25) 






3-4 


0.5807(19) 


0.5883(23) 








4-5 


0.5254(51) 










5-6 













Table 2: Average value of calculated from the starting Sierpihski carpets at different 
renormalization levels. It is obtained by considering 7 to 11 odd couplings at the sim- 
ulation temperature T sim = 1.4795. At the highest renormalization level, these average 
values are disregarded (denoted by the symbol " — ") because the eigen-exponent is not 
stable with respect to the number of the odd couplings. 



level 


SC(3,1,8) 


SC(3,1,7) 


SC(3,1,6) 


SC(3,1,5) 


SC(3,1,4) 


0-1 


1.80753(8) 


1.80762(3) 


1.80799(4) 


1.80861(5) 


1.81060(8) 


1-2 


1.81602(2) 


1.81621(3) 


1.81673(3) 


1.81857(4) 




2-3 


1.81713(1) 


1.81759(3) 


1.81948(4) 






3-4 


1.81807(1) 


1.81985(2) 








4-5 


1.81967(3) 










5-6 













16 



no 



(11 



12 



(13 













1-) 




• 


• 


2 0.5000 
































2-) 




• 




2 0.3536 








• 
























3.) 


• 




• 


2 0.2500 
























• 








4-) 






• 


4 0.2236 
























• 








5.) 


• 






1 0.2197 
























• 




• 




6.) 








4 0.1916 
























• 








7-) 






• 


4 0.1860 
































8.) 








2 0.1768 








• 
















• 


• 


• 




9.) 


• 






8 0.1734 













• 


• 




• 




• 














• 








• 


• 




• 












• 




• 




• 




• 





• 


• 






• 








• 



8 0.1655 



4 0.1613 
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